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The statistical theory of polymers tethered around the inner surface of a cylindrical channel has tradition¬ 
ally employed the assumption that the equilibrium density of the polymers is independent of the azimuthal 
coordinate. However, simulations have shown that this rotational symmetry can be broken when there are 
attractive interactions between the polymers. We investigate the phases that emerge in these circumstances, 
and we quantify the effect of the symmetry assumption on the phase behavior of the system. In the absence of 
this assumption, one can observe large differences in the equilibrium densities between the rotationally sym¬ 
metric case and the non-rotationally symmetric case. A simple analytical model is developed that illustrates 
the driving thermodynamic forces responsible for this symmetry breaking. Our results have implications for 
the current understanding of the polymer behavior in cylindrical nanopores. 


I. INTRODUCTION 

Nanoscale pores are abundant in biology and are 
rapidly gaining importance in technological and biomed¬ 
ical applications.^ Artificial solid-state nanopores enable 
the detection of single molecules via changes in ionic cur¬ 
rents through^or across^ them. Nanopores can also 
be used as molecular filtration devices, with applications 
from separating biomolecules to water purification.^ Ad¬ 
ditionally, nanoporous devices have been used as novel 
drug delivery devices^. 

Transport selectivity is one of the key challenges in 
the application of nanopores. While size exclusion is rea¬ 
sonably straightforward, chemical selectivity is harder to 
achieve. This is exemplified in nanopore-based DNA se¬ 
quencing. Recent technological advances have improved 
sensing resolution down to the single nucleotide level, 
but significant challenges remain. In particular, high fi¬ 
delity and chemically selective sensing depend on an ac¬ 
curate control of the speed with which macromolecules 
translocate through the nanopore.Interactions be¬ 
tween the DNA and the nanopore surface can lead to 
significant variation in translocation times,though 
well defined distributions have been recorded for shorter 
DNA fragments. 

While there is scope for improvement using natural 
ion channels^^ or by inserting natural protein pores in 
artificial nanopore devices,a more generic method of 
nanopore functionalization is through the coating of its 
surface with specific macromolecules. For example, bet¬ 
ter control of single-molecular transport may be obtained 
by functionalization of the pores with one-end grafted 
polymers.Solid-state nanopores grafted with poly¬ 
mers are also suggested as a possible flow control mecha¬ 
nism in microfluidic systems, responding to stimuli such 
as temperature or 

A fascinating example of a polymer-coated nanopore 
can be found in the living cell, where the nuclear 
pore complex (NPC) is responsible for all macro- 


molecular transport between the cell nucleus and the 
cytoplasm.The walls of the NPC are lined with 
natively unfolded proteins, which are intrinsically dis¬ 
ordered and rich in the phenylalanine-glycine (FG) re¬ 
peat amino-acid sequence. Macromolecular transport is 
mediated by untethered globular macromolecules known 
as nuclear transport receptors via a mechanism that is 
still disputed, though clearly dependent on the affinity of 
these receptors to the FG-repeats. NPG transport prop¬ 
erties have also been mimicked in artificial devices by 
coating solid-state nanopores with selected proteins from 
the NPG central channel. 

Besides their technological and biological importance, 
polymer-grafted pores are also of great fundamental in¬ 
terest. In the confinement of a cylindrical nanopore, 
polymer/polymer interactions have been shown, via nu¬ 
merical simulations,to yield a wide and rich pattern 
of possible morphologies. 

The understanding of the possible conformations of 
polymers within such geometries would be an important 
advance in the theory of polymer-coated pores and could 
guide experimental intuition. In particular, whilst work 
has been done to address the conformations of polymers 
within nanopores, they almost all use the simplifying as¬ 
sumption that the density of the polymers is rotationally 
symmetric around the central axis of the nanopore. In 
this paper, we address the nature of the polymer confor¬ 
mations if this restriction is lifted. 

In general, tethered polymer systems have differing 
phase behavior depending on the intermolecular inter¬ 
actions. These phases are categorized by the scaling be¬ 
havior of the radius of the gyration of the polymers^^, 
and labeled as “mushrooms” or “brushes”. When the 
polymers are tethered to the inside of a curved surface, 
such as the circumference of a nanopore, these behav¬ 
iors are complemented further by two gross phases that 
differ from the simple swelling and collapsing of a poly¬ 
mer brush. They can broadly be categorized as “wall 
phases” where the majority of the polymer density is 
found closer to their tethering points, and a “central 
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phase” where the polymers stretch away, at an entropic 
cost, from their tethering points and meet in the central 
channel of the nanopore^^’^^. Such behavior could pro¬ 
vide a gating mechanism for macromolecular transport 
through nanopores, which might be exploited in artifi¬ 
cial pores, and explain^^ the remarkable abilities of the 
NPC to selectively transport cargos as large as viruses. 

The primary driver of this change in conformation is 
the minimization of the surface area of the polymer con¬ 
densate. However, the wall phase can further reduce its 
surface area by condensing in separate clumps along the 
circumference of the pore, such that the average position 
of the polymers is still close to the wall, but they form 
into separate regions of large polymer density arranged 
around the nanopore circumference. Such clumping will 
have an impact upon the phase diagram of polymer con¬ 
formations by lowering the free energy of the wall phase. 

In order to understand this better, we determine how 
rotational symmetry breaking affects the stability of the 
previously observed central phase^^’^^. We do this using 
several methods: Monte Carlo simulations, density func¬ 
tional theory calculations and a simple analytical treat¬ 
ment. 

II. METHODS 

Our study builds upon on a previous paper^^ where we 
observed rotational symmetry breaking in Monte Carlo 
simulations of a polymer-coated cylinder. For simplicity 
and computational tractability we restrict our analysis 
to two dimensional systems, i.e., here we consider only 
variations as a function of radial and azimuthal coordi¬ 
nates. The model system is as follows: a number M of 
polymers consisting of N identical disks tethered around 
the circumference of a circle of radius i?, as schematically 
illustrated in Fig. lA. The disks are separated from each 
other by a bond length 6, and the polymers interact with 
each other through the following potential: 

00 |r|<d 

-eexp(-(|r| - ci)/(T) |r|>(i 

where d is the diameter of the disk, r is the vector con¬ 
necting the centers of both beads, e is the potential depth 
and (T represents the range of the interaction. The system 
was studied using Monte Carlo, density functional the¬ 
ory and analytical calculations. The Monte Carlo simula¬ 
tions were performed using the Metropolis algorithm^^ on 
systems of tethered polymers and the density functional 
theory is adapted from our previous work^^. The analyt¬ 
ical work is based on simple representations of polymer 
behavior. 

Monte Carlo simulations of 2D polymers with a fixed 
bond length constraint have difficulties with convergence. 
If we consider three consecutive disks along a polymer 
there are only two positions that the central disk of the 
three can take that preserve the fixed bond length con¬ 
straint. Because there is only one possible move per disk 




FIG. 1. A) We consider M polymers of N disks tethered on 
the inner surface of a circle, illustrated here for the case M — 2 
and iV = 4. B) Basis of the analytical model used to interpret 
the results. We divide the pore into Nc sections, in each of the 
tethered polymers stretch to form a clump in the centre of the 
section, as indicated by the green arrows. In this case Ac = 4. 
The change in total free energy is given by the entropic cost 
of extending the polymers together with the total energy of 
the polymer phase which depends on the change in the length 
of the interface it forms with the surrounding medium. 


for any conformation that keeps the bond length con¬ 
stant, Monte Carlo performed in this way is slow to con¬ 
verge. The disk at the end of the polymer is the only one 
which has a continuous range of possible moves. Moves 
that rotate a greater number of disks around a fixed disk 
could be considered, but the difficulty of avoiding disk 
overlap limits the likelihood that updates which move 
a large number of disks would be accepted. Therefore 
the Monte Carlo results presented here should be viewed 
more as indicative of the equilibrium structure rather 
than exact quantitative representations of it. 

A more robust but necessarily simplified density func¬ 
tional theory of the system was also used to obtain the 
equilibrium densities of the polymers. This involves con¬ 
structing a free energy functional of the polymer density 
and finding the structure that minimizes it. We adapt 
the functional used previously^^ to the here studied 2D 
system; more details can be found in appendix A. 

By considering the change in entropy of the polymers 
as they stretch away from their tethering points together 
with the the total energy of each clump, we can write 
down a simple form for the change in free energy upon 
breaking the the rotational symmetry of the wall phase. 
The total entropy change of the system will be repre¬ 
sented as a sum of terms quadratic in the distance be¬ 
tween each polymer tether point and the center of its 
clump. We can evaluate the energy of each clump by in¬ 
tegrating the potential Eq. (1) over the total area of the 
clump. This allows us to write the form of the change in 
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free energy against clump number: 

/27r(i?2 + (i?-i?j2) 
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For a more accurate investigation, we adapt our previ¬ 
ously published density functional theory^^ to describe 
two-dimensional polymers within a circular pore. Den¬ 
sity functional theory can typically provide results with 
less computational effort than Monte Carlo. In addi¬ 
tion, free energies can be obtained directly. We perform 
calculations with and without the rotational symmetry 
constraint. By comparing the resulting phase diagrams, 
the effect of rotational symmetry breaking on the phase 
of the system can be determined. We verify density func- 


where Nc is the number of clumps, Rc is the radial dis¬ 
tance the polymers stretch away from the inner surface of 
the pore, which is related to the radius of a single clump 
Rci (see Fig. IB). The only free parameter in this sim¬ 
ple analytical model is the polymer disk density p, and 
since this parameter varies when we modify e and a, we 
take this parameter from the theoretical results obtained 
in density functional theory. Next, the clump number in 
equilibrium is deteremined by minimizing AF versus Nc. 
Further details about the derivation and assumptions of 
the model can be found in appendix B. 


III. RESULTS AND DISCUSSION 

The calculations were performed with the parameters 
M = 25, = 100, b = 1.0 nm, d = 0.25 nm and R = 2b 

nm. The behavior of the system was established for a 
range of the interaction parameters e and a. 

Monte Carlo simulations can illustrate the propensity 
of the system to display rotationahy asymmetric behav¬ 
ior. Fig. 2 gives an example of the clumping behavior 
of polymers in a pore.^^ The system breaks rotational 




FIG. 2. Monte Carlo simulations of polymers tethered on a 
circular ring on the inside of a 3D cylinder (A) and to the in¬ 
side of a circle in 2D (B). The Monte Carlo simulations show 
that clumping behavior can emerge in both two and three di¬ 
mensions for interaction parameters (A) e = 0.1kbT,a = 1.00 
nm (B)e = 0.04A:6T, cr = 1.25nm, motivating more quantita¬ 
tive methods for determining the clump phase. 


symmetry for certain parameters in both two and three 
dimensions. 

Convergence issues with Monte Carlo, discussed in sec¬ 
tion III, potentially limit the reliability of this approach. 



FIG. 3. Equilibrium densities of polymer disks for density 
functional theory and Monte Carlo, for e = O.Olk^T and a = 
0.75 nm (A and B respectively) and for e = 0.04k6T and 
a = 1.25 nm (C and D respectively). Both methods yield 
similar results. 


tional theory results for the predicted clump number by 
Monte Carlo simulations for the same parameters. 

Because of the problems in the convergence of Monte 
Carlo, precise quantitative agreement between density 
functional theory and Monte Carlo is hard to achieve. 
Nevertheless, there is a qualitative similarity in the ob¬ 
tained results from the two different methods, as can be 
seen from Fig. 3. 

Fig. 4 shows an example of the various phases obtained 
from the density functional theory approach. It is ob¬ 
served that there is a multiplicity of wall states, with dif¬ 
ferent numbers of clumps. There exist many metastable 
states for each parameter set, necessitating the calcula¬ 
tions to be repeated many times to determine the phase 
with the lowest free energy. The central and the wall 
phase can both exist at the same parameter set. 

The wall phase can consist of any number of clumps 
up to the number of polymers in the system. Despite 
the differing numbers of clumps, it appears that the po¬ 
sition of the clumps reflect maximal symmetry, in the 
sense that when there are six clumps they appear to be 
positioned on the vertices of a hexagon; when there are 
five clumps they lie on the vertices of a pentagon and so 
forth. Additionally, all the clumps are similar in terms 
of disk density and spatial extent. 












4 



FIG. 4. Density in a circle of radius i? = 25 of 2D constrained polymers with parameters M = 25, b = Inm, N = 100, d = 0.25nm 
and various e and a (A) e = O.OSkbT^a = 0.75 nm, (B) e = O.Od/c^T, a = 0.75 nm, (C) e = 0.05A:6T, a = 0.75 nm, (D) 
e = 0.06kbT,a = 1.25 nm, (E) e = 0.07/c6F, cr = 1.75 nm, (F) e = O.ObkbT^a = 1.25 nm, obtained from density functional 
theory. There exist both the central phase (F) and various wall phases (A-E) with different numbers of clumps. 


According to the analytical model described in ap¬ 
pendix B, the total entropy cost is, very approximately, 
inversely proportional to the clump number. Therefore, 
large numbers of clumps carry little entropic cost of for¬ 
mation compared to a single central clump. Again ap¬ 
proximately, the total energy cost increases logarithmi¬ 
cally with clump number (up to an asymptote). The 
clumping number is determined by the delicate balance 
between these energetic and entropic contributions to the 
free energy. 

Using the density functional theory method we can test 
whether the rotational symmetry assumption has a large 
impact on the threshold interaction parameters at which 
the system has to choose between central and wall states. 
We can then investigate how well the simple analytical 
model can describe what is occuring in the phase transi¬ 
tion. 

Fig. 5 shows the phase diagram of the polymers in 
terms of central/wall phases for the cases with the ro¬ 
tational symmetry assumption (A) and without it (B). 
We can see clearly that relaxing this assumption leads to 
a wall phase that persists over a larger range a and e. 
However, the overall effect is not too large, indeed a shift 
of around 0.25 nm in a would lead to phase diagrams 
that are very similar. The analytical model (C) captures 
the behavior, despite its many simplifications. 

The success of the analytical model confirms that the 


clumping behavior results from a competition between 
the entropy of stretching the polymers against the energy 
penalty of creating a larger surface area. When there are 
large numbers of clumps, the polymers do not need to 
stretch very far, as the clump is likely to be close to the 
polymer tethering point, and therefore the entropic cost 
of forming more numerous clumps is lower. However, a 
larger number of clumps will have a greater surface area 
compared to fewer clumps, which will increase the energy 
penalty. This competition determines the equilibrium 
phase, and as we decrease the strength of the attraction 
the system will begin to favor states with more clumps. 
The analytical model gives a simple way of quantifying, 
for a particular set of parameters, how this behaviour 
emerges. 


IV. CONCLUSION 

We have shown that a system of polymers tethered 
around the circumference of a nanopore can break rota¬ 
tional symmetry, leading to the formation of clumps at 
the wall. The number of clumps can be predicted from 
the system parameters, such as the interaction strength 
and the radius of the pore. Previous theoretical results 
showed that tethered polymers could adopt wall or cen¬ 
tral phases, and this phase description remains valid even 
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FIG. 5. Phase diagram for 2D polymers with A) imposed rotational symmetry assumption, B) broken rotational symmetry 
and C) the analytical model. The blue regions are phases where the majority of the polymer density lies near the wall and the 
red regions are where the majority of the polymer density can be found near the central region. 


if the rotational symmetry assumption is dropped. When 
the rotational symmetry assumption is relaxed, the wall 
phase divides into multiple clumps which leads to a slight 
lowering of the free energy of the wall phase compared to 
the symmetric case, shrinking the boundaries of central 
phase stability on the phase diagram. 

This effect will be largely dependent on the parameters 
chosen. For instance, cylindrical pores of small radii with 
long polymers will likely always be completely clogged 
by polymer, regardless of rotational symmetry breaking. 
This is similar to what is to be expected in three di¬ 
mensions, such that previous results obtained in three 
dimensions are not qualitatively altered. 

We have developed a simple analytic model of the sys¬ 
tem that compares favorably with density functional the¬ 
ory calculations. By considering only the simplest form 
contributions to the free energy might take, we derived 
an equation that gives the change in free energy as the 
number of clumps is varied. This simple form of the free 
energy can serve to guide intution. 
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(3F[W{r)] = -^lnZ[VF(r),ro] - f VF(r)p(r) dr 
ro 

+ J y(r)p(r)dr + j dr {$"^(r) + $‘=''(r)} 
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where F is the total free energy of the system. W{r) is 
a “mean-field” potential, which can be mapped to a disk 
density p(r). I/(r) is the external potential that acts on 
the polymers, here set to 0 . gives the excess free 

energy density arising from the excluded volume imposed 
by the hard disks, and is given by^"^ 

$^^(r) = p(r)(^^^^-ln(l-i?(r))) , (A2) 

where 77 (r) is the packing fraction in two dimensions, 
given by p(r) 7 rd ^/4 with d the diameter of the disk. 
is the “chain connectivity” term, describing excess free 
energy due to the fact the hard disks are joined together. 
It is given by^^’^^ 
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Appendix A: Density Functional Theory 


Eor density functional theory calculations we use the 
following free energy functional, slightly modified from 


4>^^(r) 


1-N 

N 


In 


( 2 - ri{r) \ 


(A3) 


We introduce a Green’s function G(ro,r, 5 ; [W]) that 
satisfies 


dG{ro,r,s; [W]) 
ds 



W{r)^ G(ro,r,s; [W^]) 

(A4) 


where b is the segment length of the polymer. The 
Green’s function describes the probability that a polymer 
of s segments will extend from fq to r in the presence of 
the potential W. This equation is to be solved for every 















tethering point Tq of the polymers around the circumfer¬ 
ence of the circle. From G(ro,r,5; [W]) we calculate the 
density in the following way: 
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( r,N-s; [W])G{r, r', s; [W]) 

/dr'E.oG(ro,r',7V;[W^]) 

(A5) 

thereby taking into account polymers attached to every 
tethering point tq in the system. Z\W{y)^yq] represents 
the partition function of a polymer tethered at yq and 
interacting with the mean field. It is given by: 

Z[W{r),ro] = J G{ro,r,N-,[W])dr . (A6) 

The method for obtaining equilibrium densities re¬ 
quires the minimization of the free energy F[W{y)]. We 
use a steepest descent scheme involving the functional 
derivative of F: 


dW{y) _ 6 F[W{y)] 
dt ~ SW{y) 


the center of a clump from the wall of the pore as 
Rc, where Rc depends on Nc and is given by: 


Rc 



Nc = l 
Nc>l 


(Bl) 


• The system maintains “maximal symmetry” for in¬ 
stance, if there are 3 clumps, a rotation of the sys¬ 
tem by ^ will leave the system invariant. In other 
words, the clumps will he on the vertices of a reg¬ 
ular polygon. For practical purpose we divide the 
circle into Nc equal sectors, and all polymers teth¬ 
ered within one section form one clump at the cen¬ 
ter of that section. 

• The density of disks within a clump is constant, 
and is related to the inverse of the disk diameter 
squared. 

• The entropy difference between states is given by 
the specific entropic costs of extending the poly¬ 
mers from their tether points to the center of each 
clump. 


where we have introduced a fictional time variable t. Dis¬ 
cretizing this equation leaves us with the following form, 
which is useful for calculations: 

H.-"+‘(y) = ir"(y) + A(SM (A8) 

The complete method for solving for equilibrium den¬ 
sities involves the following: the mean field is initialized 
to a certain profile IFo, from which Eq. (A4) and Eq. 
(A5) are used to obtain the density p. Once the density 
is obtained, the functional derivative of the free energy 
at this mean field and density is calculated, and this is 
then used to update the mean field according to equa¬ 
tion (A8). This process is repeated until the system has 
reached convergence, i.e. the free energy and density no 
longer change after updates. This process is repeated ten 
times for each parameter set in order to find the global 
minimum for each case. It is important to note that the 
field Wq is initialised to random values using a random 
number generator. This ensures that the rotational sym¬ 
metry is broken: if this were not imposed, the algorithm 
would maintain the symmetry, as in our previous study^^. 


Appendix B: Analytical Model 

We denote the number of clumps as Nc and we make 
the following assumptions about them and the rest of the 
system: 

• The clumps are circular with a radius Rd. If the 
clumps lie against the wall their centers will be a 
distance Rd from the wall. If there is only one 
clump it will lie in the center of the circle. The 
clumps are identical. We represent the distance of 


• The entropic contribution to the free energy of 
clump formation from a uniform wall configuration 
will be represented by a sum of terms of the form: 


A5i(r) 



(B2) 


for each polymer where r is the end to end distance, 
kb is the Boltzmann constant. 


The system being treated is shown in Eig. 1. In this 
case, Nc = 4:. The polymers stretch from each quadrant 
edge towards the center of that quadrant’s clump. 

The free energy change of this system is given by: 


AF = AE- TAS (B3) 


Where T is the temperature of the polymer. The total 
entropy change AS will be given by the entropic cost of 
extending the polymers: 


TAS = -Nc^ ^(i?^ + (i?-i?e)'-2(i?-i?c)i?cos(^,)) • 

i=l 

(B4) 

This can be understood as a summation of the contribu¬ 
tions of all the polymers in a sector as they stretch from 
r = R and 0 = Oi to r = R — Rc and 0 = 0, where 9 is 
the angle along the sector. This is illustrated in Eig. 6. 
Eor analytical reasons, we assume that the polymers are 
tethered uniformly and continuously along the circum¬ 
ference of the circle, in which case we would modify our 
discrete representation to: 


TAS{Nc) = -Nc 


hT 



2{R-Rc)R cos{e)) de 
(B5) 
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The total internal energy of a clump is given by the 
sum of interdisk interactions within the clump. We as¬ 
sume that the inter-clump interactions are negligible. 
This assumption will hold when there are few clumps, 
such that the distance between clumps is large. 

-Be = ^ / /3(r')/3(r)9^(i' “ r')drdr' (BIO) 

To calculate the internal energy, we integrate over the 
potential in a clump using a radial coordinate system. 
The full form of the inter disk interaction is: 


(j){Y — r') = eexp 




e(lr-r'l-d) 


(Bll) 


where 0 is the Heaviside step function and |r—r'| is given 
by (in plane polar coordinates): 

|r —r'l = (r^ -h r'^ — 2rr' cos{q — q'))^ (B12) 


FIG. 6. The geometry used to calculate the entropy of ex¬ 
tension of the polymers. The polymers tethered around the 
edge of the pore stretch to the center of the clump, at ^ = 0 
and r = R — Rd- In order to calculate the entropy change, 
contributions proportional to the square of this extension are 
summed over polymer tether points lying on the edge of the 
circle. 


where ps is the number of polymers tethered per unit 
length of pore circumference, i.e 27Tps = M 
This integral can be evaluated to yield: 


TAS{Nc) 


^ 271 Nh‘^ 

(27i{R^^{R-R^f) 


^AR{Rc — R) sin 



For example when A^c = 1 we use the definition of Rc 
given before to give the total entropy change as: 


TAS{Nc = 1) = 


MhTR^ 

Nb‘^ 


(B7) 


We use the incompressibility condition to calculate the 
clump radius: 


such that 


NM _ At] 
NcTtRIi ~ ^ 


Rci 



(B8) 


(B9) 


Where 77 is the packing fraction, as defined in density 
functional theory. 


where r and r' are the radial positions, relative to the 
centre of the clump, of points r and r' and q and q' are 
the angular coordinates with respect to the clump center. 
Using the assumption that the disk density in the clump 
is constant and rearranging we obtain: 


/ d\ r^ci r‘^7T 

Ec = ep^TT exp ( — ) / dr dr' dqrr' exp 

V^/Jo JQ Jo 


X 0(|r-r'|-^^) 


(B13) 


where we have performed the integration over one of the 
angles. If we assume that the diameter of the bead is 
small in comparison to the radius of the clump (again 
true for low numbers of clumps) then we ignore the Heav¬ 
iside function for analytical tractability, giving: 


Ec = TTcp^ exp ( — I Ref dx fdx' [ dqxx' exp 

J Jo Jo Jo 



(B14) 

where x = r/Rc and x' = r'jRc. The integral can be 
evaluated approximately to give: 


Ec = ep^exp ^ —^ ^1 —exp 

(B15) 

and the total energy will be given by = NcEc. Using 
the relation between Rc and Nc allows us to write the 
energy in terms of clump number. 


E =eexp 



NMcF^TTp X 


(B16) 


— exp 


( 1 I nm\ 
[~a\l^p) 




















such that the energy change as a function of the number 
of clumps is given by: 


AE" = —eexp 



NMcF^TTp X 


(B17) 


exp 


( 1 I nm\ 




The combined entropy and energy are given in Eq. (2) 
in the main text. An example of the free energy change 
against clump number is shown in Fig. 7. 


AF [A:,T] 



FIG. 7. Form of the analytic free energy change against the 
clump number Nc for a specific set of parameters: N = 
100, M = 25, d = 0.25 nm, R = 25 nm, b = 1 nm, e = 
0.07 k^T, a = 0.5 nm with these parameter sets the packing 
fraction p is 0.383989. 


Using similar arguments, we are able to derive a form 
of the free energy for the equivalent system in three di¬ 
mensions: 


= 


3MkbT + 4^ sin 




(B18) 


l2r]MNa^ee‘^/‘^ (2 - + 2 ^) 

with Rc = when Nc > 2 and Rc = R when 

Nc = l. 
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